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Abstract. 

In the present work we discuss based on extensive fully atomistic molecular 
dynamics simulations the dynamics of Neon atoms oscillating inside (5,5)singlc- 
walled carbon nanotubes (CNTs) and boron-nitride (BNNT) ones. Our results show 
that sustained high frequency oscillatory regimes are possible for a large range of 
temperatures. Our results also show that the general features of the oscillations are 
quite similar to CNT and BNNT, in contrast with some speculations in previous work 
literature about the importance of broken symmetry and chirality exhibited by BNNTs. 



PACS numbers: 68.65. k, 61.46. +w, 68.37.Lp, 71.15. m 
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1. Introduction 

Since the pioneering work of Cumings and Zetttl pQ carbon nanotube (CNT) based 
oscillators have became object of many studies. Cumings and Zettl demonstrated the 
controlled and reversible telescopic extension of multiwalled CNTs, thus realizing ultra 
low-friction nanoscale linear bearings. Zheng and collaborators [21 E] have proposed that 
this phenomenon could be exploited to built nano-oscillators in the frequency range of 
gigahertz. Legoas et al. [H Ej carried out molecular dynamics (MD) simulations and 
concluded that sustained oscillations are possible when the radii difference between 
inner and outer oscillator moving parts are about ~ 3.4 A. Frequencies as high as 
40 GHz are possible and even clear chaotic signatures were observed [6]. Rivera et 
al. [7] studied the behavior of oscillator damping via MD and they showed that the 
damping depends on the inverse of the CNT length, which is in good agreement with 
the prediction of analytical models [21 [3]. They also showed [8] that the frequencies in 
the gigahertz decrease as the length of the tubes increases. Zhao et al. [H] studied the 
oscillator energy dissipation energy mechanism, also using MD simulations, and they 
demonstrated the importance of the overlap length of double CNTs and the mechanical 
deformation of the outmost CNT to this mechanism. In spite of the large amount 
of investigation of CNT based oscillators some aspects remain not clearly understood, 
especially with relation to friction mechanisms [101 EJ EE] • In this sense it is important 
to study other related systems such as oscillators made of different materials and also 
atoms and molecules oscillating inside the tubes. Only very few studies have been 
reported along these lines, such as oscillators composed of boron nitride nanotubes 
(BNNTs) [13], BNNTs in association with CNT ones [Hj, C 60 molecules [US M\ and 
rare-gas atoms p21 [18] inside CNTs. The choice for using rare-gas atoms is associated 
with their inert nature which prevents chemical reactions with the tubes. It has been 
demonstrated [13, [18] that sustained oscillations are possible for a large range of energy 
variations and also, in some circunstantes, the tube chirality is important to determine 
the oscillatory patterns. However, there are conflicting data for the behavior of the 
oscillator frequency decaying with the tube diameter. Further studies are necessary to 
better understand these aspects. In the present work we report extensive fully atomistic 
molecular dynamics simulations for the study of Ne atoms oscillating inside (5,5) CNTs 
and (5,5)BNNTs, to our knowledge this is the first study of rare-gas atoms oscillating 
inside BNNTs. 

2. Methodology 

We have carried out molecular dynamics simulations using the universal force field UFF 
[19] [20] as implemented in the Cerius2 package |21j. This force field includes van der 
Waals, bond stretch, bond angle bend, inversion, torsion and rotation terms. It has been 
used with success in the study of dynamical properties of carbon based nanostructures 
[H El [221 1231 121] • Initially, we perform a molecular mechanics energy minimization to 



Neon Atoms Oscillating Inside Carbon and Boron-Nitride Nanotubes 



3 



optimize the geometry of each isolated nanotube. Then, a Ne atom is placed about 
lA from the nanotube ends and the dynamical simulations are carried out. We have 
used in all microcanonical ensemble (NVE). We have used an integration step 

of 1 fs. The MD simulations were performed in the following temperatures (in Kelvins): 
5, 10, 20, 40, 80, 160, 240, and 300. In order to thermically equilibrate the system a 
running of 200 ps is initially performed. After the equilibration a running of 50 ps is 
carried out and the data trajectory (positions and velocities) saved. From the trajectory 
data we perform a fast Fourier transform(FFT) to obtain the oscillatory frequencies. We 
have considered CNT and BNNT with (5,5) chirality type. Two different lengths were 
considered (25 A and 50 A in order to investigate the frequency dependence with the 
length of the tubes. The choice of (5,5) nanotubes is due to that they have the minimum 
diameter that allows the Ne atom encapsulation and also this particular chirality implies 
fewer atoms in the unit cell which reduces substantially the computational cost of the 
simulations. 

3. Results and Discussions 

In table [T] we show the calculated oscillation frequencies for a Ne atom inside a CNT, 
for temperatures from 5K up to 240K. For 300K we observed that the Ne atom is not 
encapsulated. In our simulations all the atoms of the nanotube are free to vibrate and 
with the increase of temperature the amplitude of the vibrations of the atoms at the 
end of tube increase. Mainly due to the van der Walls interactions a dynamical barrier 
is created which prevents the Ne atom encapsulation. As we can see from the Table 
stable frequencies are possible for a large range of temperatures, which could be an 
interesting aspect to be exploited in possible applications, such as clocks. For the 25 
A CNT from 5 up to 240 K the frequencies are 136.71 and then they drop to 117.18 
GHz. For the longer tube (50 A), as expected, this occurs at lower temperature (160 
K), the frequencies drops from 78.12 to 58.59 GHz. This sustained oscillatory regime is 
better visualized in Figs. 1 and 2. These results are consistent with the ones obtained 
by Zeng et al. [H] in spite of their different protocol simulations. They used an impulse 
dynamics that can be mapped into our NVE simulations varying the initial conditions 
with respect to the temperature of simulations. 



Table 1. Calculated frequencies (GHz) for Ne@CNT(5,5). 



Temperature (K) 


CNT 


5 


10 


20 


40 


80 


160 


240 


25 A 


136.71 


136.71 


136.71 


136.71 


136.71 


136.71 


117.18 


50 A 


78.12 


78.12 


78.12 


78.12 


78.12 


58.59 


58.59 



In Table [2] we present the equivalent results for the case of a Ne atom inside a 
(5,5)BNNT. Again, for the cases investigated, the temperatures where the Ne atom is 
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Figure 1. Molecular dynamics simulation results for Ne atom oscillations inside a 
CNT of 25A length for the following temperatures a) 80K, b) 160K and c) 240K 

not encapsulated were 300K for the tube of 25 A, and 240 and 300 K for the tube of 
50 A. We believe that this occurs by the same reasons (large end atom tube movements 
creating a dynamic barrier) as in CNT case. That this occurs at a lower temperature 
in the BNNT case can be attributed in part to different symmetries (buckling) and the 
slightly diameter differences between the CNT (6.8 A) and BNNT (6.9 A). We expect 
larger fluctuations for tubes of larger lengths due to mass differences. 



Table 2. Calculated frequencies (GHz) for Ne@BNNT(5,5). 



Temperature (K) 


BNNT 


5 


10 


20 


40 


80 


160 


240 


25 A 


136.71 


136.71 


136.71 


136.71 


136.71 


117.18 


117.18 


50 A 


78.12 


78.12 


78.12 


78.12 


78.12 


78.12 
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Figure 2. Molecular dynamics simulation results for Ne atom oscillations inside a 
CNT of 50A length for the following temperatures a)80K, b)160K and c)240K 

Again, as observed for the CNT case, sustained oscillatory regimes are possible for 
a large range of temperatures, and at almost the same frequency values. The only major 
differences are at what temperature we observe dropping frequency occurring values and 
the absence of encapsulation at 240 K for the BNNT of 50 A. This could be in part 
attribute to mass and diameter tube differences. The major point however it is that 
in contrast to some previous works [T7] for the systems investigated here the symmetry 
breaking and chirality of BNNT do not seem to play important roles in determining the 
oscillatory regimes. The general oscillatory behavior (Figs. 3 and 4), frequency values 
and axial displacement movements (Fig. 5) are quite similar to the ones observed in the 
CNT case. 

In summary, we have carried out detailed analysis of extensive fully atomistic 
molecular dynamics simulation for Ne atom oscillating inside carbon and boron nitride 
nanotubes of (5,5) chirality and different lengths (25 and 50 A). Our results show that 
well defined sustained oscillatory regimes are possible for a large range of temperatures 
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Figure 3. Molecular dynamics simulation results for Ne atom oscillations inside a 
BNNT of 251 length for the following temperatures a)80K, b)160K and c)240K 

(from 5 up to 240 K). There are several stable frequencies for each one nano-oscillator 
configuration studied here, ranging from 58 to 137 GHz. Another interesting aspect of 
these nanostructures is that we can tune the frequency by a suitable modification of the 
nanotube length and/or type. Our results showed that the results for the BNNT are 
quite similar to the CNT ones in terms of general oscillatory behavior, frequency values, 
and axial displacement movements. In this sense our results are closer to the obtained 
by Zeng et al. [18j and in disagreement with the ones obtained by Lin et al. [17] . We 
hope the present work to stimulate other studies along these lines. 
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Figure 4. Molecular dynamics simulation results for Ne atom oscillations inside a 
BNNT of 50A length for the following temperatures a)40K, b)80K and c)160K 
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